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ABSTRACT 

We report the first detection of C-C3H2 in a circumstellar disk. The C-C3H2 
J=6 — 5 line (217.882 GHz) is detected and imaged through Atacama Large Mil- 
limeter Array (ALMA) Science Verification observations toward the disk around 
the Herbig Ae star HD 163296 at 0.8" resolution. The emission is consistent 
with that arising from a Keplerian rotating disk. Two additional C-C3H2 tran- 
sitions are also tentatively detected, bolstering the identification of this species, 
but with insufficient signal-to-noise ratio to constrain the spatial distribution. 
Using a previously developed model for the physical structure of this disk, we fit 
a radial power-law distribution model to the C-C3H2 6 — 5 emission and find that 
C-C3H2 is present in a ring structure from an inner radius of about 30 AU to an 
outer radius of about 165 AU. The column density is estimated to be 10^^-10^'^ 
cm~^. The clear detection and intriguing ring structure suggest that C-C3H2 has 
the potential to become a useful probe of radiation penetration in disks. 

Subject headings: protoplanetary disks; astrochemistry; stars: formation; ISM: 
molecules; techniques: high angular resolution; radio lines: ISM 
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Introduction 



Pre-main sequence stars are commonly surrounded by disks, which serve the dual pur- 
pose of funneling accretion onto the central star and transporting away angular momen- 
tum. These disks are also the formation sites of planets and connect protostars and planets 
physically and chemically. The main component of gas-rich circumstellar disks, cold H2, 
cannot be observed directly because of a lack of transitions at suitable energy levels. In- 
stead, constraints on disk mass, density and temperature structures depend on observations 
of trace species, mo st commonly dust and CO gas (e.g. iQi et al.l l2006l : iPietu et al.l 120071 : 
Andrews et al.ll2009l ). To access other disk properties, such as ionization levels. X-ray and 
UV radiation fields, turbulent mixing, chemical age and bulk composition requires the de- 
velopment of additional molecular probes whose formation and destruction depend critically 
on one or several of these properties. 

The past decade has witnessed significant progress on disk chemistry theory, resulting in 
large numbers of predictio ns on how molecular abundance structures should depend on dif- 



ferent disk 


properties (e.g. 


Ai 


cawa 


Foeel et al. 


2011 




Walsh et al. 


2012) 



^awa k Nomurall2006l: IWillacvl 120071 : ISemenov fc Wiebell2011 



2OI2I ). A lack of observational sensitivity has however limited 



detections to the simplest molecules in the millimeter /submillimeter range which character- 



^^^^ ^^^^ ^^^^ O"" 

CS, so, CHo, N2H+ and HC.N ( 


3utrev et al. 19971: Aikawa et al. 


2OO3I: ' 


Thi et al. 


2004: 


Semenov et al. 20( 


5: Dutrev et al. 


2007: Oi et al. 


20081: Fuente et al 


2010: 


Hennine et al. 


2010: Obere et al. 


2010. 2011). This list is expected to expand with the arrival of ALMA 



and its unprecedented sensitivity and spatial resolution, which will enable the development 
of a more diverse set of molecular probes. In this Letter, we use ALMA Science Verification 
observations of HD 163296 to report the first detection of cycloprop enyhdene, C-C3H2, in 
a circumstellar disk. At a distance of ~122 pc (jPerryman et al.l 119971 ). this Herbig Ae star 
(stellar mass 2.3 M0; spectral type Al; age ~4 Myr) harbors a large disk (a scattered light 



pattern extending out to a radius of ^500 AU. iGradv et a. 



tinuum and molecular l i ne emission ([Mannings fc Sargent 



2004: 


Isella et al. 


2007: 


Oi et al. 


2011) 



1997 



20001) with strong millimeter con 



Thi et al. 2004: Natta et al. 



make this system an ideal target to search for new molecules. 

The c-C.sH-;. ni o lecule is a small cyclic hydrocarbon, which was first detected in space 



by iThaddeus et al.l (Il985l ) towards Sgr B2. Since then, C-C3H2 has been detected in a wide 



tarv nebulae and extraealactic sources fe.e.. Matthews & Irvine 1985: Cox et a . 


1988, 1989: 


Madden et al. 1989: Benson et al. 1998: 


Menten et al. 


1999: Fosse et al. 2001: 


Sakai et al. 


2008: Tenenbaum et al. 


2009: Liszt et al. 


2OI2I), and is 


one of the most abundant molecules 
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with 3 carbon atoms in the interstellar medium (ITeyssier et al.l |2004| ) . In the context of 
disks, C-C3H2 has some favorable characteristics that may make it a powerful molecular 
probe, especially when observed in combination with isotopologues and chemically related 
molecules: 



C-C3H2 has two equivalent H nuclei, which couple to generate ortho (nuclear spin of 1) 
and para (nuclear spin of 0) states. Deviations from the statistical weights of 3 to 1 
may be used to constrain the formation conditions of the molecule (jPark et al.l 120061 ) . 



C-C3H2 is also one of two isomers; the other is the linear molec ule /-C3H2. The is omer 
ratio can be used to assess the fractional ionization in clouds (IFosse et al.l 1200 ll ) and 
may provide similar constraints on disks. 

The formation and destruction chemistry of C-C3H2 (and carbon chains) are sensitive 
to high energy radiation. The vertical and radial abundance structur e should provide 
stron g constraints on the penetration depth of X-rays and UV photons (ISemenov fc Wiebe 



2OIII). 



Within the family of hydrocarbons, the sensitivity to turbulent mixing increases dra- 
matically from C2 to C3H2, and from C3H2 to larger hydro carbons. Observed ratio s 
and limits have the potential to constrain transport in disks (ISemenov fc Wiebell201ll ). 



Finally, the chemistry that produces C-C3H2 is important to constrain in its own right, 
since C-C3H2 shares precursors with other, more complex carbon chains, which may be an 
important source of carbon during planet formation. In light of these possibilities, we aim 
to provide first constraints on the abundance, distribution and o/p ratio of C-C3H2 in the 
HD 163296 disk. 



2. Observations 

HD 163296 (RA: 17^56""21!287, DEC: -2r57'22'.'39; J2000.0) was observed on 2012 
June 9, June 23 and July 7 as part of the ALMA Science Verification (SV) program in band 
6. At the time of observations, 25 12-m antennae were used (2 7-m ACA antennae were 
present in the dataset but flagged) in the array with a resulting synthesized beam size of 
0.9x0.7" (PA=84°). The source was observed with a total integration time of 84 minutes. 
The ALMA correlator was configured to simultaneously observe four spectral windows (two 
in each sideband). One of those windows (SpwID ^1) covered both ^'^CO and C^^O 2 — 1 lines 
(220.399 and 219.560 GHz, respectively). The data were calibrated in the CASA software 
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package (v3.4), following the detailed calibration and imaging scripts provided by the ALMA 
science verification team. Since those scripts, along with fully calibrated measurement sets, 
are publicly available onlin^, we do not repeat the details here. The 218 GHz (1.37mm) 
continuum was generated from the line-free channels in the two spectral windows (SpWID 
7^0 and #1) of the lower sideband, and the fiux is determined to be 608.5±2.5 mJy, which 
agrees with the SMA observations (IQi et al.ll201ll ). Excellent agreement is also seen between 
the SMA and ALMA spectra for the ^^CO and C^^O 2-1 lines with fluxes within 10% 
Joi et allbonl ). 



All of the C-C3H2 lines are continuum-subtracted and retrieved from SpwID #0 with a 
channel width of 0.488 MHz (0.67 km s~^) covering the total 1.875 GHz bandwidth from 
216.167 to 218.042 GHz. We have binned the lines in channels of 1 km s~^ to achieve higher 
signal-to-noise ratios and still keep the key kinematic features of the disk. The resulting line 
rms noise level is estimated to be 2.2 mJy beam^^ per 1 km s~^ channel in line free channels. 



3. Results and Analysis 

3.1. C-C3H2 detection 

Figure [1] shows that the C-C3H2 6 — 5 line is clearly detected towards HD 163296 and 
that the position angle of its resolved velocity fleld is consistent with that of ^^CO from the 
same data set, which shows a clear pattern of Keplerian rotation. Interestingly, the C-C3H2 
6 — 5 moment map displays a double-peaked morphology, indicative of an inner hole. The 
C-C3H2 3 — 2 and 5 — 4 lines are only tentatively detected, insufficient to resolve the velocity 
field or any spatial structure. To improve the sensitivity in imaging these lines, we applied 
a Gaussian taper to the visibility data that corresponds to 1" FWHM in the image domain, 
which resulted in a beam of 1.3x1.2" (PA=78°). 

Because the line emission changes position with channels, using a single mask for all 
channels to extract spectra will lead to a degradation of the signal. Instead, we create 
separate masks to cover the C-C3H2 6 — 5 line emission at each channel and use these masks 
to extract the line spectra. The same masks are also used for other C-C3H2 lines. Figure 
|2] shows the resulting spectra of all three C-C3H2 3 — 2 lines together with ^'^CO 2 — 1 line 
(with its own masks covering the ^'^CO 2 — 1 emission at each channel). The 6 — 5 line 
is clearly detected with the expected line shape and Visr consistent with ^^CO. The 5 — 4 
line presents a consistent spectral profile, at lower signal-to-noise ratio. In contrast, the 



^https:/ / alniascience.nrao.edu/alniadata/sciver/HD163296Band6 
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3 — 2 line is both narrower and clearly blue-shifted, suggestive of a different physical origin, 
perhaps in a disk wind. However, both the 3 — 2 and 5 — 4 lines are observed at the edge 
of the spectral window, so this apparent asymmetry should be confirmed by independent 
observations before any strong conclusions are drawn. 

The integrated fluxes are reported in Table [H These are just below the upper limits 



repor ted in two previous failed searches of C-C3H2 in disks ( iFuente et al.ll2010l : lOberg et al. 
2010h . 



3.2. Ring structure model 



We explore the distributions of C-C3H2 in HD 163296 based on a previ ously develope d 
accretion disk model with well-defined temperature and density structures (IQi et al.ll201ll ). 
To simulate and compare with the data, we assume the disk material orbits the central 
star in Keplerian motion, and fix the disk geometric and kinematic parameters that affect 
the observed spatio-kinematic behavior of the disk. The details of the stellar and accretion 
properties and the disk geometric and kinematic parameters are summarized in Table 3 of 
Qi et al.l ( 120111 ). We have developed a physically self-consistent accretion disk model with 
an exponentially tapered edge that matches the spectral energy distribution and spatially 
resolved millimeter dust continuum emission. The disk temperature and density structures 
are further constrained by multiple CO and CO isotopologue line observations. Such analysis 
provides the essential framework for constraining the distribution of molecular abundances 
in protoplanetary disks. 

Within this model framework, we assume that C-C3H2 is present with a constant abun- 
dance in a layer with boundaries toward the midplane and toward the surface of the disk 
( Oi et al.ll2008l: lOberg et al.ll2012l ). This assumption is motivated by chemical models (e.g. 



Aikawa &: Nomura 20061 ) that predict a three-layered structure where most molecules are 



photodissociated in the surface layer, frozen out in the midplane, and have an abundance 
that peaks at intermediate disk heights. The surface (ag) and midplane (0"^) boundaries 
are presented in terms of S21 = T,h/{^-59 x lO^^cm"^), where T^h is the hydrogen column 
density measured from the disk surface in the adopted physical model. This simple model 
approach approximates the vertical location where C-C3H2 is most abundant. Due to the lim- 
ited signal-to-noise of the multiple C-C3H2 line detections, we cannot constrain the location 
of this vertical layer based on the C-C3H2 data. We therefore fix the vertical surface bound- 
ary ds as 0.32 (logio(o"s) = —0.5) and the midplane boundary am as 3.2 (logio(o"m) = 0.5), 
between which lies th e expected location of na ost molecules with peak abundances in the 
warm molecular layer ( lAikawa fc Nomurall2006l ). 
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We model the radial column density distribution of C-C3H2 as a power law Nioo x 
[r/lOOy, where Nioo is the column density at 100 AU in cm~^, r is the distance from the 
star in AU, and p is the power-law index. Since the emission shows a clear indication of 
a central hole, we fit for both an inner radius Kin and outer radius Kout together with 
the power-law parameters (Nioo and p). Using the structure model and the fixed vertical 
distributions, we compute a grid of synthetic C-C3H2 J=6 — 5 visibihty datasets over a range 
of Rout, Rin, P and Nioo values and compare with the observations. The best-fit model is 
obtained by minimizing x^, the weighted difference between the real and imaginary part of 
the complex visibility measured in the {u, f )-plane sam pled by the ALMA observations. We 



use the two-dimensional Monte Carlo model RATRAN (iHogerheijde fc van der Takll2000l ) to 



calculate the radiative transfer and molecular excitation. The collisional rates are taken from 



Molecular Database ( jSchoier et al.ll2005l ) 



Chandra &: Kegell (I20001) and the raolecu lar data file is retrieved from the Leiden Atomic and 



Figure |3] shows the surfaces for the Ri„ and Rout versus the power law index p. We 
find the p is constrained between -2.5 and -1.5, R^ 15-40 AU and Rout 150-200 AU (within 
Icr) with best-fit values at -2, 30 and 165 AU, respectively (reduced = 1-07). By fixing 
the above parameters at their best-fit values, N(c-C3H2) is determined to be 2.2±0.2 xlO^^ 
cm~^ at 100 AU. Figure H] presents comparisons between the observed channel maps and the 
best-fit model. The residual image doesn't show any significant emission, which indicates 
that the C-C3H2 emission is consistent with that arising from a Keplerian rotating disk. 
The 6 — 5 line is actually blended with the 61^6 ~ 5o,5 (ortho) and 60,6 — 5i^5 (para) lines at 
the same frequency. Si nce the Einstein A- coefficients for both ortho and para lines are the 



same, 5.393x10"^ s~^ (jSchoier et al.l 120051 ). we fit the data with the ortho line only and the 
resulting column density can be accounted as the sum of both ortho and para forms of the 
molecule. 

For the weaker 5 — 4 ortho line, the signal-to-noise ratio is too low to provide any 
constraints on the C-C3H2 the distribution. By assuming the C-C3H2 distribution derived 
from the 6 — 5 line, however, we can use the unresolved line flux to constrain the o/p ratio. 
The resulting best-fit N(C3H2-(o))=1.6±0.3 xlO^^ ^nd N(C3H2-(o + p)=2.2±0.2 xlO^^ ^^^-2 
at 100 AU (as shown above). Hence o/p=2.7±1.7, which is close to the statistical value of 
3, but the large uncertainty imphes that better data are required to use this ratio as a probe 
of formation conditions. 
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Discussion 



We have clearly detected C-C3H2 in the disk around HD 163296. This is one of the 
two largest molecules detected in disks so far, the other being HC^iN, also a carbon -chain 
molecule, detected recently using deep single-dish observations (IChapillon et al.ll2012l ). 



Using the spatially and spectrally resolved line emission, we have constrained the radial 
distribution of C-C3H2 in the HD 163296 disk as a ring structure ranging from ~30 to ~165 
AU. The outer edge of <250 AU (2a limit) is much smaller than the gas disk size which 
extends to 500 AU (IQi et al.ll201ll ). which suggests that C-C3H2 formation is slow or that 
C-C3H2 destruction is fast in the outer disk. In terms of formation, the best-fit outer radius 
coincides with the previously determined location of the CO midplane snow line, and it is 
possible that CO freeze-out limits the carbon available in the gas phase to form hydrocarbons 
in general and C-C3H2 in particular. A second characteristic of the outer disk is that it is 
rathe r tenuous because of a rapidly decreasing column density outside of 200 AU ( lOi et al. 
2OIII ). This may result in efficient penetration of radiation and thus destruction of the easily 
dissociated C-C3H2. Explicit model predictions are required to resolve which effect drives the 
disappearance of C-C3H2 in the outer disk. 

No existing disk chemistry model in the literature contains predictions for the radial 
distribution of C-C3H2. Based on strong correlations f ound between C-C3H2 and CCH i n 
several astrophysical environments ( iLucas fc Lisztll2000l : iPety et al.ll2005l : iGerin et al.ll201ll ). 
predictions on CCH in disks can be used to shed further light on the origins of C-C3H2. 
Aikawa &: Herbst (j2001 ) predicts a CCH central cavity in the absence of X-ray radiation ; 



when X-rays are present , and HD 163296 is a known X-ray emitter (jSwartz et al.l 12005 



Giinther fc Schmittlboogf ). the CCH column is instead centrally peaked. Despite the presence 
of X-rays, the central cavity could be due to a UV dominated radiation field. This may be 
tested by observing C-C3H2 towards T Tauri stars with weak UV fields and strong X-rays, 
and towards additional Herbig Ae stars with weak X-rays. 

The inner cavity could also be a product of the main formation pathway of C-C3H2. Car- 
bon chains can form through at least three different pathways: (1) ion-neutral reactions at low 
temperatures, (2) photo-erosion of larger carbonaceous compounds, and (3) ne utral- neutral 
react i ons at luke-wa.rm temperatures following CH4 ice evapo ration at ~30 K (IHerbst et al. 
19841 : iTevssier et aPbooi ISakai et al.lboosi ICerin et al.lboill ). Detailed modehng of C-C3H2 
radial distributions in each of these formation scenarios are clearly needed to use C-C3H2 
rings as tracers of the radiation field and other disk characteristics. 



In addition to more detailed modeling, additional observations are needed to constrain 
whether different hydrocarbons, especially CCH and C-C3H2, are correlated in disk envi- 
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ronments, and to constrain the vertical distribution of these molecules. The latter requires 
multiple lines with a range of excitation energies. With the combined advance in theory 
and observations, C-C3H2 may become one of the more useful probes of penetration of the 
radiation fields in disks. 

Facility: ALMA 

This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00010.SV. 
ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS 
(Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with 
the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and 
NAOJ. 

We would like to thank an anonymous referee for thoughtful suggestions for the paper. We 
also acknowledge NASA Origins of Solar Systems grant No. NNX11AK63. 
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Table 1: C-C3H2 line results. 



Transition 


Frequency (GHz) 


(K) 


Beam 


J Fdv (niJy km s ^) 


33,0 ~ 22,1 


216.279 


19 


I'.'S X l'.'2(76°) 


53 [9] 


61,6 ~ 50,5/60,6 — 5i,5 


217.822 


39 


0'.'9 X 0'.'7(83°) 


185[10] 


5l,4 - 42,3 


217.940 


35 


1'.'3 X l'.'2(78°) 


74 [9] 
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Fig. 1. — The integrated intensity maps summed between and 11 km s~^) and intensity- 
weighted mean velocity fields of ^^CO 2 — 1 and C-C3H2 6 — 5 lines (left panel), C-C3H2 5 — 4 
and 3 — 2 lines (right panel) toward HD 163296. The resolved velocity field of the C-C3H2 
6 — 5 line agrees with the CO kinematics. In the C-C3H2 maps, the first contour marks 3a 
followed by Icr contour increases. The rms varies between 6 and 9 mJy km s~^ per beam. 
Synthesized beams are presented in the lower left corners. The star symbol indicates the 
continuum (stellar) position. The axes are offsets from the pointing center in arcseconds. 
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V,sr [km s-^] 

Fig. 2. — Extracted spectra of the three C-C3H2 transitions toward HD 163296, plotted 
together with ^^CO 2 — 1. The dashed hne marks Vlsr toward HD 163296. 
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Fig. 3. — Iso-x^ surfaces of i^out and Ri^ versus p. Contours correspond to the 1-5 a errors. 
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Fig. 4. — The top panel shows the velocity channel map of the C-C3H2 6 — 5 emission 
toward HD 163296 (velocities binned in 1 km s~^). Contours are 0.0022 Jy Beam^^ (la) 
X [3, 4, 5, 6, 7, 8, 9]. The middle and bottom panels show the best-fit model and the difference 
between the best-fit model and data on the same contour scale. The star symbol indicates 
the continuum (stellar) position. The axes are offsets from the pointing center in arcseconds. 



